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Abstract 

In this paper we present for the first time a complete description 
of the Bohm model of the Dirac particle. This result demonstrates 
again that the common perception that it is not possible to construct 
a fully relativistic version of the Bohm approach is incorrect. We ob- 
tain the fully relativistic version by using an approach based on Clif- 
ford algebras outlined in two earlier papers by Hiley and by Hiley and 
Callaghan. The relativistic model is different from the one originally 
proposed by Bohm and Hiley and by Doran and Lasenby. We obtain 
exact expressions for the Bohm energy-momentum density, a relativis- 
tic quantum Hamilton- Jacobi for the conservation of energy which in- 
cludes an expression for the quantum potential and a relativistic time 
development equation for the spin vectors of the particle. We then 
show that these reduce to the corresponding non-relativistic expres- 
sions for the Pauli particle which have already been derived by Bohm, 
Schiller and Tiomno and in more general form by Hiley and Callaghan. 
In contrast to the original presentations, there is no need to appeal to 
classical mechanics at any stage of the development of the formalism. 
All the results for the Dirac, Pauli and Schrodinger cases are shown to 
emerge respectively from the hierarchy of Clifford algebras C13, C30, Cqi 
taken over the reals as Hestenes has already argued. Thus quantum 
mechanics is emerging from one mathematical structure with no need 
to appeal to an external Hilbert space with wave functions. 



'E-mail address b.hiley@bbk.ac.uk. 
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1 Introduction. 



In a recent paper, Hiley and Callaghan flj have shown that both the Schrodinger 
and Pauli theories can be completely described within a pair of nested Clif- 
ford algebras, Co,i and £3,0 and there is no need to use the wave function. 
Thus contrary to popular belief there is no need to regard the Hilbert space 
as an indispensable basic feature of quantum mechanics. Instead, the infor- 
mation required to reproduce all quantum effects can be carried within the 
algebra itself, specifically by appropriate elements of a minimum left ideal 
as first pointed out by Frescura and Hiley [2] where the reasons for adopting 
this approach were discussed. 

We showed in our previous paper PQ that the Clifford algebraic formal- 
ism was completely equivalent to the conventional approach to quantum 
mechanics. This opens up the possibility of a different interpretation, an 
interpretation that brings us closer to providing an explanation of quantum 
phenomena in terms of a non-commutative geometry. Furthermore, and 
somewhat surprisingly, we found our formulation enabled us to make con- 
tact with the Bohm approach, casting a completely different light on this 
interpretation. Specifically in the case of the Schrodinger theory, we ob- 
tained exactly the formalism used in the Bohm model [3], [I] and [5]. While 
for a Pauli particle, which of course carries spin, we found a generalisation 
of the Bohm, Schiller and Tiomno [6], [7] theory. Both these models are 
non-relativistic and therefore the previous paper contained no new physics, 
being of a pedagogical nature and written to show explicitly the role played 
by the Clifford algebra. The purpose of this paper is to apply the same 
methods to the Dirac theory. 

In conventional relativistic quantum mechanics, the Clifford algebra made 
its appearance indirectly as an attempt to remove the negative energy that 
arises in the relativistic expression for the energy, E = ±^p 2 + m 2 . (We use 
natural units throughout). Dirac proposed that we start from the positive 
linear Hamiltonian H = a.p + (3m where a and j3 were unknown elements 
that have to be determined by the constraints placed on the energy. Thus in 
order to satisfy E 2 = p 2 + m 2 , a and (3 must be four anti-commuting objects 
that could be represented by 4 x 4 matrices. With the linear Hamiltonian, 
the equation of motion becomes 



where $ is a four-dimensional column vector, assumed to be an element in 
a finite dimensional Hilbert space. This can be simplified by writing the 



idt$> = [a.p + f3m}^ 
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equation in manifestly covariant form 

(i'fdp - m)tf = (1.2) 

This reinforced the belief that the wave function was the natural vehicle for 
describing quantum phenomena. However what Dirac had actually discov- 
ered was that a and f3 were elements of the Clifford algebra Ci^. The Pauli 
spin matrices had already alerted us to a possible role for the Clifford al- 
gebra, but the Schrodinger theory seemed not to require a Clifford algebra. 
However once one realises that Cq i = C then we see that the Schrodinger 
theory can also be discussed in terms of a Clifford algebra. The way the 
Clifford structure can be used in the Schrodinger and Pauli approaches to 
quantum theory was described in detail in a previous paper pQ, The pur- 
pose of this paper is to show how these Clifford algebraic techniques can be 
applied to the Dirac theory. In order to bring out the quantum aspects of 
our approach, we will confine our discussion to the free, charge-neutral Dirac 
particle. It is easy to extend our approach to charged particles by simply 
introducing minimal coupling to the electo-magnetic field in the usual way. 

2 Previous Approaches. 

We will assume that the application of the conventional approach through 
spinors represented in a Hilbert space is by now well known. However at- 
tempts to discuss the Dirac theory from within the Clifford algebra itself 
have been judged to have achieved limited success. Here we are thinking 
specifically of the work of Eddington [8J, Takabayasi [9] and Hestenes [TO] 
and references therein. In this paper we will extend our previous work [T] 
to the Dirac theory. We will also show that, as a consequence, we are able 
to provide a complete relativistic generalisation of the Bohm approach. 

Previous attempts to build a Bohm model have generally been restricted 
to using the expression for the Dirac current to obtain trajectories for the 
Dirac electron. This approach has been used by Bohm and Hiley [4], |12j and 
Holland [5], where some applications have been discussed. Gull, Lasenby and 
Doran [13] have also used the Dirac current to calculate trajectories in their 
investigation of quantum tunnelling in the relativistic domain. However in 
no case was a quantum Hamilton-Jacobi equation for the conservation of 
energy derived and consequently no relativistic expression for the quantum 
potential has been found. A very early proposal was made by Bohm him- 
self [TT], but this proposal was found to be unsatisfactory as will become 
apparent from the work presented in this paper. These partial attempts 
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have been seen as a critical weakness of the Bohm model and has led to 
the general conclusion that the Bohm approach cannot be applied in the 
relativistic domain. This paper shows adds further evidence to show that 
this conclusion is not correct. 

In view of the arguments proposed in [I] and in [2J, we show how we 
can also dispense with the wave function \E' in the Dirac theory and instead 
use an element of a minimal left ideal to carry the information normally 
carried by the wave function. Thus once again we can work entirely within 
the algebra with no need to introduce an external Hilbert space structure. 
This, of course, not only provides us with an alternative approach to the 
Dirac theory itself, but it also provides a way to generalise the Bohm model 
so that it can be applied to all relativistic particle situations. This work 
complements the work of Horton, Dewdney and Nesteruk |15j |16| |17| who 
discuss the application of the Bohm approach to relativistic scalar particles 
and vector bosons. 

To obtain our results we must examine the Clifford algebra in more 
detail and abstract all the information necessary to describe a relativistic 
quantum process. We should point out that previous attempts at exploiting 
this Clifford structure have been made by Hestenes [18], [ID]. The method 
we use in this paper is more direct, being based on the ideas discussed in 
Hiley [19] and in Hiley and Callaghan pQ. 

3 The Clifford Algebra C 1)3 . 

3.1 C13 over the reals. 

In order to treat the Dirac theory in the same way as we have treated the 
Schrodinger and Pauli theories, we need first to understand some of the 
complexities introduced by the larger Clifford algebra C\^. This algebra is 
generated by {1,7 M } where \i = 0,1,2,3 and [7^,7^] = 2g tiu . We choose the 
diagonal metric (1,-1,-1,-1). In our approach we are taking the Clifford 
algebra as basic and projecting out a space-time manifold using the mapping 
77 : 7^ — > e M . Here e M is a set of orthonormal unit vectors in a vector space 
V\ 3 , the Minkowski space-time for an equivalent class of Lorentz observers 

3.2 The Primitive Idempotents 

The next step is to examine the primitive idempotent structure in £1,3, as 
these play a key role in our approach. In the case of the Schrodinger and 
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Pauli algebras, the primitve idempotents have a very simple structure. For 
example in the Schrodinger algebra, the primitive idempotent is trivial, be- 
ing the identity of the algebra. In the Pauli algebra all primitive idempotents 
can be generated from e = (1 + crj)/2 where Oi is one of the generators of 
the algebra C%$. These idempotents form an equivalence class under the 
transformation e' = geg" 1 where g is an element of the spin group SU(2). 

In the Clifford algebra C\ 3, taken over the reals, we can distinguish three 
classes of primitive idempotents [20j . characterised by 

ei = (l + 7o)/2; e 2 = (l+ 7 3o)/2; e 3 = (1 + 7m)/2. (3.1) 

Each of these idempotents can generate an equivalence class using the rela- 
tions e' { = ge%g~ x where, in this case, the g are elements of the spin group 
SL(2C). Crumeyrolle [20J calls the elements in each set 'geometrically equiv- 
alent', since they are equivalent under the Lorentz group. The three defin- 
ing idempotents can be related to each other through where 
Uij(^ g) is some other invertible element of the algebra. 

The reason for choosing one specific idempotent from the set (|3.ip is 
determined by the physics we are interested in. For example, the choice 
of (1 + 7o)/2 is made in order to pick out a particular Lorentz frame by 
fixing the time axis defined by 70 through the projection rj defined in a 
previous section. Any other Lorentz frame can be obtained by choosing 
the appropriate Lorentz transformation, g = A(v) from Sl(2C) and forming 
7q = A _1 (v)7oA(v). In this way we relate the algebra to an equivalence 
class of Lorentz observers. 

Alternatively (l+73o)/2 can be chosen if we want to pick out a particular 
spin direction, the 3-axis in the Lorentz frame defined by 70. In this case 
we are highlighting the even sub-algebra C° 3 = €3^. Other spin directions 
can be described by transforming 730 using appropriate elements of the spin 
group. 

If we want to consider the little group of the Lorentz group, 50(2, 1) 
[21], then either (1 + 7123) /2 or (1 + 7 3 o)/2 can be chosen. These elements 
project onto the idempotents of the corresponding Clifford algebra £2,1- 

As we have remarked above, the three primitive idempotents themselves 
are equivalent under a larger group, the Clifford group which consists of 
all invertible elements of the algebra. For example we find the following 
relations 

e 2 = (l-73)ei(l + 73)/2 
e 3 = (1 - 75M1 + 7s)/2 
£3 = (1 - 7012)62(1 + 7oi2)/2 
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This group contains reflections and inversions and is thus clearly related to 
parity and time reversal. We will not discuss further these relations in this 
paper. 

3.3 Formation of elements of a minimal left ideal in C lj3 over 
the Reals. 

We start by choosing the primitive idempotent e\ = (1 + 7o)/2. Thus we are 
working in a particular Lorentz frame. A general element of the minimal 
left ideal is generated by forming A(l + 7o)/2, A being all the elements of 
C\ t 3. This minimal left element can then be written in the form 

2$ L (^) = o(a^)(l + 70) + 6(a^)(7i2 + 7012) + c(^)( 723 + 7023) 
+d(x^)(7i 3 + 7013) + /(x M )(7oi - 7i) + 2(2^X702 - 72) 

+M^)(703 - 73) + n(^)( 75 - 7123) (3.2) 

Here the set {a(x^), b(x^), c(x^), d(x^), f{x^), g(x^), h(x^), n(x^)} are the 
eight real functions that can be used to specify the quantum state of the 
Dirac particle. Recall that this element encodes all the information that 
is normally encoded in the wave function which uses 4 complex numbers. 
Thus there exists a relation between these two sets of parameters which will 
be given in section 3.7, equation (|3.1ip . 

Examining (13.2p shows that it splits into an even parlQ, § e L , and odd part, 
with the same set of parameters occurring in each expression. Since the 
same information about the system is contained in each part, it would be 
convenient if we could use only one of these expressions in describing the 
state of the process. We will follow Hestenes |22| and use only & e L which, 
for convenience, we will write as (pL- Specifically we write 

4> L = a + 6712 + C723 + c?7i3 + /7oi + 5702 + ^703 + "-75 (3.3) 

The set {a,b,c,d,f,g,h,n} is determined by solving the Dirac equation as 
was first shown by Sauter [23] . 

We can put these relations in a simplified and more transparent form by 
writing the basis of this ideal as 

[1,712,723,713,701,702,703,75] (l + 7o)- 
Then we can split this expression into two parts 

[1,712,723,713] (l + 7o) (3.4) 
by 'even' we mean that <& e L only contains even products of the 7s. 
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and 



[75,701,702,703] (l + 7o) = 75 [1,712, 723, 713] (l + 7o) (3.5) 

It should be noted that the square bracket in f)3.4[) and the second square 
bracket in (|3. 5j) are the same and clearly generate the quarternion algebra. 
Thus the elements of the ideal split into two quaternions, explaining the old 
name 'bi-quternion' used by Eddington [8j. Splitting the element in this 
fashion immediately enables us to distinguish between its large and small 
components. To see this recall that an infinitesimal Lorentz boost, Sv, along, 
say the x-axis, can always be written in the form A(5v) = (1 + jqi5v) so that 
when 5v — > the contribution of the 701 is negligible. Thus, in general, when 
the boost is negligible, the contribution of the 70^ terms will be negligible, 
hence the term 'small'. One we can show that this division corresponds 
exactly to the conventional division of the Dirac wave function into its large 
and small components. 

3.4 Real or Complex field? 

There is a difficulty in taking over the reals which we must now ad- 
dress. In the Schrodinger and Pauli cases we noted that there was no need 
to introduce the complex numbers in order to compare with the standard 
approach. In both these cases, the Schrodinger Clifford, Co,i, and the Pauli 
Clifford, C3 0, contain an element in the centre of the algebra that squares to 
minus one. In the case of Co,i it was, trivially, the generator of the algebra 
e, while for the Pauli, £3,0, it was the element ei23, the product of all three 
generators. The Clifford algebra £1,3, only has the identity in the centre, 
therefore we cannot simply use C\ 3 over the reals. This would make it im- 
possible to reproduce all the results of the standard theory which exploits 
C®Ci >3 . 

However if we go to the next larger Clifford algebra, which is C^\- We find 
this algebra does contain a non-trivial term in the centre, namely, -F01234 = 
i*5, where Fa {A = 0, 1 . . . 4) are the generators of the algebra £4,1, a notation 
introduced by Rocha and Vaz [23]. We will now show how we can still use 
the real number field and yet account for all the properties of the standard 
approach. To do this we will exploit the isomorphism = C ® Ci 3. This 
will enable us to use taken over the reals. The structure is now rich 
enough for us to be able to relate the eight real parameters {a, . . . , n} to the 
four complex numbers, 1^1,1^2^3 and 1/J4, contained in the standard Dirac 
vector in Hilbert space (i.e. the Dirac spinor). 
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In C^i the generators Fa satisfy the anti-commutaion relations [Fa,Fb]+ = 
IGab where Gab is the metric tensor with diagonal elements (—1, 1,1, 1,1). 
Here the element F§ = .Foi234) being in the centre of the algebra, will play the 
role of the imaginary i. Specifically we use the projection p : — > <C®C\^ 
defined by, 

p(i ? M)=7M ^ = 0,1,2,3. (3.6) 
p(F 5 ) = i (3.7) 

One significant point to notice is the way the signature changes from (4, 1) 
to (1,3). This is characteristic of the way the next lower Clifford algebra 
is embedded in the even part of the higher algebra. This result follows 
immediately from the nature of the projection p as we will now show. 
Consider the anti-commutator 

[F^4, F U 4] + = —2G^ U p = 0,1,2,3. 

Then forming 

[p(F^),p(F M )] + = -2p(G^) 
and using equations (|3.6p and (|3.7p we find 

= -2G Mi/ = 2g [lv 

where is the metric of C\^. This shows why the the signs of the metric 
alternate. 

Let us now consider the Pauli algebra which is an even sub-algebra 
of Ci 5 3 under the projection 

A(7fco) = °k 
^(75) = C123 

In this case C^ 3 = giving the metric tensor as g-ij = (1, 1, 1). Following 
this through to the lowest algebra Co,i we have the sequence of embeddings 

£2,4 — ^4,1 —> Ci,3 £3,0 — > Co, 2 — > Co,i- 

Here we have chosen to start this sequence with the conformal Clifford C2A 
since this contains the Penrose twistors [25]. Physically this means we have 
a hierarchy of Clifford algebras which fit naturally the physical sequence: 

Twistors —> relativistic particle with spin — > non-relativistic particle with 
spin — > non-relativistic particle without spin. 
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3.5 The Dirac algebra 



Let us now examine the structure of C 4j i in more detail. Again we find it 
contains three classes of primitive idempotent characterised by (1 + .F 4 )(1 + 
*bl4)/4, (l + i ? 04)(l + i ? oi4)/4 and (l + i ? i 4 )(l + i ? i234)/4. All other primitive 
idempotents can be derived from these using elements of the spin group. Out 
of all the possible idempotents we choose the specific idempotent 

e F = (1 + F 04 )(l - F 34)/4 = (1 + F 04 )(l + F 5 F 12 )/4 

we obtain the 
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If we now use the projection p defined in (I3.6P and (J3T 
primitive idempotent that we will be working with, namelj 



e 7 = (l + 7o)(l + i7i2)/4. (3.8) 

The reason for our choice now becomes apparent. The first term in e 7 allows 
us to project out a shadow manifold, which as we have shown, is equivalent 
to picking out a particular Lorentz frame by choosing a given time axis. 
The second term allows us to pick out a spin plane which is perpendicular 
to the spin axis 63 . This will enable us to see exactly how the Pauli and the 
Schrodinger Clifford sub-algebras are contained in the larger algebra. This 
particular choice of idempotent allows us to use directly the relation of the 
hierarchy of Clifford algebras discussed in the previous sub-section. Having 
said this, it is important to realize that, mathematically, we could choose 
any primitive idempotent, but we have chosen (13.8P because it enables us to 
compare our results directly with those of the conventional approach using 
the Dirac representation as we will show in section 3.7. 



3.6 The Minimum Left Ideal Generated by e 7 . 

At this stage we need to relate the minimal ideal given in (13. 3p to the 
minimum left ideal generated by the idempotent e 7 (|3.8p . It can be shown 
by straight forward multiplication that the basis of the minimal left ideal 
generated in this way is 

e 7 = (1 + 7o + «7i2 + «7oi2)/4; 

723^7 = (723 + 7023 + «7l3 + «7013)/4; 

703£ 7 = (703 - 73 + 175 - «7i23)/4; 
7oi e 7 = (701 - 7i - *702 + «72)/4. 

2 It looks as if we have introduced the complex number i. However we are actually 
using F5, which for short hand convenience, we are calling i. 
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To make contact with the element of the ideal shown in equation (|3.2p . we 
write 



4$l 7 = (a- ib)e 1 + (c - ic0723e 7 + (h - m)7 30 e 7 + (/ + 10)70167. (3.9) 

To check the relation of this element with that given in (j3.2[) we must expand 
out this expression and examine the even part. We find immediately that 
$L = 25J[3>£ 7 ]. If we examine the imaginary part of we find that by 
multiplying from the left by 721, we have 5ft [^z, ] = 721^[ < £l 7 ]- Thus the 
same information is encoded four times in the minimum left ideal generated 
by e 7 given by f|3.8 j) . 

3.7 The Standard Matrix Representation 

Since this paper is mainly about the use of Clifford algebras in physics, we 
need to relate the information contained in elements of minimal left ideals 
with that contained in the wave function. For comparison we choose the 
standard representation used by Dirac |26j . In this representation we have 



Here we write 75 = 70123 • 

If we substitute these matrices into the equation (|3.9|) and compare the 
resulting matrix with the usual Hilbert space column vector \P with compo- 
nents ipi, tp2i fa, "04?) we fi n d 

ipi = a — ib; ip2 = —d — ic, if) 3 = h — in, ip4 = f + ig, (3.11) 

This will enable us the write 



2a = (^+Vi); 2i6=(^-Vi); -2cf=(V>2* + V>2); 2ic=(^-^); 
2h = (r 3 +i> 3 ); 2m = (V> 3 *-V>3); 2/ = (^ + V 4 ); %g = ^-^\). 



We will use these relations later in the paper to compare our results with 
those of the standard approach. 

3.8 Minimal Right Ideals. 

In an earlier paper |19j we argued the once we have decided on a particular 
idempotent we need to construct the Clifford density element [CDE], p c = 
^L^ij, where <&r is the conjugate to There are two ways to construct 




(3.10) 



10 



this conjugate. Firstly we can use the direct method by forming the product 
(1 + 7o)x4 so that 

2§ R = a(l + 70 ) - 6(712 + 7012) - £(723 - 7023) - d(7i3 - 7013) 

-/(701 - 7i) - 5(702 - 72) - M703 ~ 73) + "-(75 - 7123) (3.12) 

Again we see this element of the right ideal splits into an even and odd 
part, containing the same information so we will once again consider only 
the even part, *ff R which we will again write as i/jr for convenience. Then 

<j> R = a- 6712 - C723 - £^713 - /7oi - 5702 - /*703 + "-75 (3.13) 

We can obtain this result more simply by taking the Clifford anti-morphism 
called Clifford conjugatiordof ^>l- Symbolically this is written as ^>r = ^l, 
where ~ denotes the anti-morphisrrQ. 



3.9 The Relation of ^r to the Dirac Adjoint Spinor 

Now let us return to the relation between the element of the right ideal 
and the adjoint spinor as defined in standard Dirac theory. We need to 
start with the element of the left ideal defined in equation (|3.9p and find 
the conjugate minimum right ideal which also involves taking the complex 
conjugate because F5 = — F5 so that §r = &* L . Thus we have 

4<S> R = (a + ib)(j)i - (c + id)<f) 2 ~(h + in)<j) 3 - (/ - ig)fa 

where 

4>i = (1 +70 + H12 + H012); 

(f>2 = (723 + 7023 - «7l3 - ^7013); 

03 = (703 + 73 + «75 + ^7123) ; 

04 = (701 + 71 + *702 + ^72)- 
If we now use the relations (|3.1ip , we find 

4<f> R = ifrlEh. + r 2 E 2 ~ r 3 E 3 - 1>IE 4 . (3.14) 

3 Any element, C, of the algebra is, in general, a sum of its scalar, vector, bivector, 
axial vector and pseudo-scalar parts, viz, C — S + V + B + A + P. The Clifford conjugate 
is given by C = S -V - B + A + P. 

4 Throughout the rest of this paper we will write for 4> l. 
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where the {4>i} have been replaced by their matrix representation counter- 
parts, {Ei} which are the row vectors E\ = (1 0); Ei = —i(0 1 0); E3 = 
(0 10) and £4 = (0 1). This then gives the expected bilinear invariant 

$R$L = |VM 2 + IV>2| 2 " |^3| 2 - |^4| 2 

In the standard Hilbert space approach this corresponds to where ip is 
the adjoint wave function, rp = t/^To- Thus we see that we can identify the 
content of the element of the right ideal that we have constructed with the 
content of the adjoint Dirac spinor. 

4 The Clifford Density Element and Bilinear In- 
variants 

4.1 The Clifford Density Element 

Having defined an element of a minimal left ideal and its conjugate right 
element, we can now form the key quantity, the Clifford density element, 
p c = &l^r = 4>L^4>R- From the discussion in the previous section, we can 
write 

fafa = (a 2 + b 2 + c 2 + d 2 - f 2 - g 2 - h 2 - n 2 ) + 2(an -bh-cf + dg)^. 

It is notationally more convenient to write 0x in the form 

4> L = R e ^l 2 U (4.1) 

where U is a bivector such that UU =1. It is not difficult to show that 
U is a general element of the spin group, SL(2C) which is, of course, the 
covering group of 50(1,3). 

The expression (14. ID is identical to what Hestenes [10J called the 'spinor 
operator'. As far as the work we a reporting here, we do not use the term 
'spinor operator' because the notion of an operator has no meaning in the 
algebraic structure we are discussing. For us it is simply the even part of an 
element of the minimal left ideal §l = i?e 7s ^/ 2 L/"(l + 7o)/2. 

In this notation, the element of the corresponding right ideal is 

4> R = XJRe^l 2 (4.2) 

so that 

p c = = (j) L {\ + 7o)0R 

= R 2 e^{l + U-/ U) (4.3) 
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With this notation we find 

§L<i>R = R 2 cos j3 + i? 2 75 sin /3 

4.2 The Bilinear Invariants 

In order to make contact with the physical parameters that describe the 
quantum process, we construct bilinear invariants directly from the CDE, 
p c . We have already seen that we define this element through 

Now we have the freedom to choose the idempotent e, but as we pointed out 
earlier, we want to compare our results with the results obtained from the 
standard approach so we will choose 

e 7 = (1 + 70 + «7i2 + «7oi2)/4 (4.4) 

The set of 16 bilinear elements written in the conventional theory are 

n = 

n = (% 5 |^) 
jm = (<j>|Y*|¥) 

2S^ = *(%^]|*} 
> = -(^l*) 

where 7^ = 17 s 7^. 

One obvious bilinear invariant we need is the Dirac current, = (^\^\^). 
In the algebraic approach this current is given by 

J" = trW<t> L e<t> R ) 

Here for simplicity we write tr for the normal trace divided by four. The 
only non- vanishing part of the trace comes from the vector part of 4>L£<t>R-, 
that is from 4>l1$4>r- Let us write 

4>LlQ<pR = ^ A V J U . 

Then 
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A tedious but straight forward calculation shows that 

A = a 2 + b 2 + c 2 + d 2 + f 2 + g 2 + h 2 + n 2 

A 1 = 2{af -bg + cn- dg) 

A 2 = 2(ag + bf-ch- dn) 

A 3 = 2(ah + bn + cg + df) 

Then using the relations between the parameters {a, b . . . n} and the compo- 
nents of the standard Dirac wave function {ipi, . . . ^4} defined in equation 
(|3,lip . we can write the Dirac current in the standard representation, viz 

J° = |Vl | 2 + I V>2 I' + |V>3 I' + |^4 1' 

J 1 = ^4+^3+^2+^4$; 

J 2 = ityrfl - + V-3V-2 - I'M 

J 3 = ^1-^1+^1-^*2 

Thus we have shown that the Dirac current can be written more simply in 
the algebraic form 

J = 4>L1°4>R- (4.5) 

This form for the current has already been used by Hestenes [27] and 
Lounesto [28j. 

In our approach the Dirac current is calculated by simply using the 
second term in the idempotent (j4.4p . It is then useful to ask what the 
other terms produce when taken individually. We find that (fti^R produces 
a scalar and a pseudo-scalar, viz, 

(t> L <t>R = + ^75- (4.6) 

A vital property for the Dirac particle is its spin. Following Messiah [29] 
we find the total angular momentum is 

J = L + S with L = r x p, S = a/2 

from which we define the spin bivector as 2S^ V = i{^>\^~f v \^) . In the 
algebraic approach this spin bivector can be written in the form 

2 P S = 4>l1vi4>r. (4.7) 
The spatial components of the spin given by s % = S^ k with i, j, k cyclic. 
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To complete the set of 16 elements we also have the dual of the current, 
JP, sometimes referred to as the axial vector current or even the Proca 
current. Algebraically this is obtained from the expression 



We have now formed the 16 bilinear invariants that are defined by the CDE, 
p c . It should be noticed that these invariants are obtained from each term 
in the idempotent so that 

n + hj 5 = <j) L <j) R . 

J = 4>LlO<t>R- 

2S = 4>l7i24>r 

Thus we see all the information about the bilinear invariants is contained in 
the CDE, p c = $^(1 + 7q + £712 + ^7012)^' R, each element in the idempotent 
giving a bilinear invariant. This means that in our approach the quantum 
state of the process is described directly in terms of these invariants and is 
therefore directly observable unlike the wave function. To find the tensor 
components to match up with the conventional theory, we have 



n = tr[(/> L <f> R ] 
J» = tr[Y<t>LlM 



If we try to specify the quantum state of the process using these 16 invari- 
ants, all observable quantities, we will find the system over specified since 
we only need 8 real parameters {a,b,c,d, f,g,h,n} to define the quantum 
state. However, as Takabasayi [9] has already pointed out, these bilinear 
invariants are not linear independent. In fact he shows that there are only 7 
independent bilinear invariants, finding 9 subsidiary conditions relating the 
invariants. We list these here for convenience. 



J = 75^L7O120i? = 4>L73<i>R 



(4.8) 



and 




fj,, v, a, A cyclic 
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where p 2 = Q, 2 +Q 2 is a positive definite scalar. A more detailed examination 
of p shows that it corresponds to ip'ip where ip is the conventional wave 
function. Thus we can interpret it as the probability of finding the particle 
at a particular position. For a more detailed discussion of the relations 
between these invariants we refer the interested reader to Takabayasi's paper 

Since we only have seven independent bilinear invariants, something is 
missing. We need one extra bit of information to completely pin down the 
quantum state. Takabayasi suggests that we need to construct a set of 
bilinear invariants of a second kind. These are defined as follows 

^roj^ = 9 M *r* - ^>Td^. 

where T is one of the sixteen independent terms formed by products of the 
generators 7**. Here we are using a short hand notation sometimes used in 
field theory (see Raymond |30|). 

Since we need only one type-two invariant to complete our description 
of the state, we have a choice. Our choice, which is different from the one 
chosen by Takayabasi, is the energy-momentum tensor 

2iT> w = 4>^(d u ip) - (d^h'V ( 4 -9) 

where ifi is the conventional wave function spinor, while tp, is the adjoint 
spinor. As we will show in the next section, this will help us identify the 
relativistic generalisation of the Bohm energy and the Bohm momentum 
that we have used in the Schrodinger and Pauli theories, showing that these 
parameters are not arbitrary but are needed to complete the specification 
of the quantum state. 

4.3 The Energy-momentum Density. 

Our first step is to write the energy-momentum tensor (|4.9j) in algebraic 
form. This then becomes 

2iT^ = tr{rW<p L )e(k R - MPfa)]} = tr^^eP 'fa)] 

Since the only non- vanishing trace is a Clifford scalar, and since 7^ is a 
Clifford vector, we must find the Clifford vector part of 4>l^P '<t>R- We are 
again using the idempotent given in equation (|3.8p . namely 

e' = (1 + 70 + 1712 + £7oi2)/4. 
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We find the only term in that gives a Clifford vector is 7012, so that we 
need only consider 

2T^ = trtftfiP-yoiriR)]- (4-10) 

Now let us first find the energy density T 00 . In order to do this we need to 
evaluate (<^l^7oi20_r) in terms of the functions {a(x^), . . . ,n(x M )}. 

After some straight forward but tedious calculations we find 

(t>iP-yoi2<t>R = AZ(xP)7i (4.11) 
where the A^ are given by 

Al = -( a Pb + cPd+ fPg + hPn) 

A\ = -(aPg + bPf + cPh + dPn) 

A\ = (aPf -bPg-cPn + dPh) 

Al = {aPn-bPh + cP f -dPg) 

If we use the relations ()3.1ip . we can put T 00 in a more familiar form in 
terms of the wave function. Thus 

4 



where ipj are the four complex components of the standard Dirac spinor and 
we have also written ipj = Rj exp iSj with Rj and Sj real functions. 

We claim this is exactly the relativistic version of pEs where Eb is the 
Bohm energy density. To show this, we find that in the non-relativistic limit 
Y2j=iRj®tSj Yl~j=i(^j^tSj). The latter is just the expression for pEs 
found in the Pauli case. (See Hiley and Callaghan pp.) Furthermore this 
reduces to Eb = —dtS, the well known expression for the Bohm energy for 
the Schrodinger particle. Thus we identify the Bohm energy through the 
relation 

P E B = T 00 . 

Similarly we can also show that the momentum density can be written in 
the form 



3 



i E(V^Vj - ^d^j) = E R2 J VS r 
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so that we can write the Bohm momentum as 

pP k = T ko 

Then it is not difficult to show that this again reduces, in the non-relativistic 
limit, to the Bohm momentum found in the Pauli case and reduces further, if 
the spin is suppressed, to the well known Schrodinger expression Pg = VS. 
This condition is sometimes known as the guidance condition, but here we 
have no 'waves', only process, so this phrase is inappropriate in this context. 

To summarise then we have defined the Bohm energy-momentum vector 
through the relation 

2pP» = 2T"° = tr[ 7 o (0L^7oi2&i)]. (4.12) 



4.4 The Appearance of Two Currents in the Relativistic 
Theory. 

In the earlier attempts to apply the Bohm approach to the Dirac theory, 
Bohm [11], Bohm and Hiley [12] and Gull, Lasenby and Doran, [13] used 
the Dirac current to provide a means of calculating particle trajectories. 
This was done under the assumption that these would be a generalisation 
of the trajectories in the non-relativistic case calculated from the expression 
for the Bohm energy-momentum density as was investigated by Dewdney et 
al. P3]. 

The Dirac current is identified with a four momentum P^ through 

p pv = m{ijyfil)) 

or in component form 

p[E,P]=pm O f /3?/0,OW) 

However to show that this is not the Bohm energy-momentum defined by 
p p£ = recall that the Dirac equation flL2]) can be written in the form 

where jjl = i^d^. We can now form 

together with the adjoint form 
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Here the Dirac current, J M = (i/j^i/j), appears on the RHS of both the last 
two equations. Thus subtracting these two equations, we find 

Now let us compare this expression with 

This expression is clearly not the same as that given for the current JP. 
However expanding the LHS of the equation for J^, we obtain the Gordon 
decomposition 

2mJ^ = i [i>{d^) - (d^)ip] 

where 2ia^ = (7^7^ — 7^7^) . This shows that the current can be decom- 
posed into a convection part, the first two terms, and a spin part. 

If we examine the convection part of the current, although it is closer 
in form to T^°, it is still different. Thus we have two distinct currents, one 
is the convection current obtained from J M and a Bohm energy-momentum 
density current derived from T^°. Both are conserved. 

This difference should not be surprising because the conserved Dirac 
current, is the Noether current obtained from a global gauge transforma- 
tion, while the energy-momentum density is the Noether current produced 
by a space-time translation. This difference does not appear in the case of 
the Pauli and the Schrodinger particles, the two currents turn out to be the 
same. 

From the physical point of view it seems, at first sight, that having two 
conserved currents is strange. Contemporary field theory emphasizes and 
uses the Dirac current J**, since this is term that couples to the electro- 
magnetic field. The energy-momentum density is simply used to show that 
when we integrate it over a volume element, it leads to the global conserva- 
tion of energy and momentum. [See for example Schewber [31].] The local 
energy-momentum density is given no physical meaning. The main point of 
the Bohm model is to give meaning to the energy-momentum density of the 
individual particle. 

Attempts have been made to exploit the difference between these two 
currents in another context, namely, an extended spinning relativistic object. 
Takabayasi [9], [52] and Bohm and Vigier [33] have tried to understand this 
difference by considering a relativistic liquid drop. Here they assume that a 



19 



collection of particles are in some form of relative motion. This implies we 
have both a circulation of particles and a circulation of energy. We would 
expect these to be different in the relativistic domain. Thus the liquid drop 
can be described in terms of the relative motion of two centres, the centre 
of mass-energy and the centre of particle density. 

In the same way we can think of the Dirac particle as an extended process 
that must be described, in the first approximation, by a process that has 
two centres of motion which coincide in the rest frame. In the conventional 
field approach, integrating the energy-momentum density over space avoids 
raising questions about any possible internal structure of a Dirac particle. 
On present experimental evidence this seems to be justified in the case for 
leptons, although when we come to baryons, we are forced to introduce inner 
structure in terms of quarks, but the quarks themselves are again assumed 
to be point-like structures leaving open any deeper structure. At one level, 
our investigation here can be regarded as simply exploring the consequences 
of this bilocal structure. Let us now move on to consider in more detail the 
structure of the energy-momentum density and its time evolution. 



5 The Time Evolution of the Energy-Momentum 
Density. 

5.1 Symmetrized Energy Constraint 

In order to investigate time development, we need to relate neighbouring 
points on the base manifold at neighbouring times. Normally this would be 
done by using single Dirac derivative. However because we are exploiting 
the structure of a Clifford bundle, there are two types of derivatives. The 
reason for needing two derivatives was explained in Hiley [19]. These are 
written in the form 

i5 = 7 ^ and 15 = 1)^ 

These derivatives give us two time development equations, one equivalent 
to the Dirac equation which when expressed in terms of an element of a 
minimal left ideal is 

i^d^L ~ m^L = (5.1) 

The other is equivalent to the adjoint Dirac equation expressed in terms of 
an element of the corresponding right ideal, 

i&^rf + m ^R = ( 5 - 2 ) 
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The energy constraint insists that we must have (d^d 11 + m 2 )ijj = satisfied 
for both and <& R , since the quantum state is described by the CDE 
p c = $l<&r. Thus two second order derivatives are involved (d^d^&^&R 
and ^L(d fl d ll ^ R ). These are used to produce two equations by taking their 
sum and their difference. The sum gives the energy conservation equation 

(d^<5> L )<S> R + $ L (d^$ R ) + 2m 2 $ L <S> R = 0. (5.3) 

While the difference produces the following equation 

^L(d^ R ) - (d^L^R = (5.4) 

which, as we will show below, describes the time evolution of the spin and 
its components. 

5.2 The Quantum Potential 

Now we will use equation (j5.3j) to investigate energy conservation. To anal- 
yse this equation further we need to see where the Bohm energy-momentum 
as defined in equation (|4.12p fits in as was done for the Pauli particle. To 
proceed let us first introduce a more general variable defined by 

2 P P» = [(3^)70120* - <pL7oi2(d^ R )} (5.5) 

Let us also introduce a quantity 

2pW» = -^(^7012^) 

Combining these two equations, we obtain 

(9^)70120* = p[P»-W»} 
and - tanaWfl) = plP^ + W] 

which can then be written as 

-Wfa = [P^-W^ L70 i2 (5.6) 

and 

d»<t> R = 7012<MP" + W] (5-7) 
Since we are going to use equation (|5.3p . we can use equation (|5.6p . to form 

-a^L = [a^ - d^w] + [P" - w»] d^ LlQ12 
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After some algebra and finally multiplying from the right by ipR, this equa- 
tion can be written in the form 

- (<W L )0 R = P [p^ + w^w* - (p»w^ + w^p*)] 

+ [d^ - d^W] (5.8) 

This gives us the first term in equation (|5.3p . Now we must consider the 
second term in this equation. Repeating an analogous set of steps but now 
using equation (|5.7j) , we find 

+(0L7oi2^) [d^ + dnW»] (5.9) 

Substituting both these equations in equation (j5.3j) . we finally find 

P 2 + W 2 + [Jd^ - d^P^J] + [Jd^W^ + d^W^J] - m 2 = (5.10) 

Here we have used the relation 2pJ = V'L7oi2V'-R> where J is essentially the 
axial current. This term reduces to the spin of the Pauli particle in the 
non-relativistic limit. Equation f|5. lOj) can be further simplified by splitting 
it into its Clifford scalar and pseudoscalar parts. The scalar part is 

P 2 + W 2 + [Jd^W^ + d^W^J] — m 2 = (5.11) 

This is to be compared with the energy equation 

Pt_tP^ — m 2 = 

Thus we see that the extra two terms must be related to the quantum 
potential in some way. Before we arrive at an exact expression for the 
quantum potential, we must first note that the momentum, P^, as defined 
in equation (|5.5p is not yet the Bohm momentum defined in equation (|4.12jl . 
Equation (I4.12D tells us that Pg is the 7 coefficient in the expression for 
P 1 *. However it is not difficult to abstract the Bohm momentum from the 
P 2 term in equation (jS.lip . To do this we need to recall equation (\5.5h and 
use equation (|4.12p to find 

3 

Using the definition of Pg given in equation (|4.12p , we find 

3 

Ap 2 P 2 = Ap 2 P 2 +Y J A ^ A i 

i=l 
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If we write 

3 

i=i 

we then find equation (|5,lip can be written in the form 

Pj + n 2 + W 2 + [Jd^W 1 * + dpW^J] - m 2 = 

Then we see that the quantum potential for the Dirac particle is 

Q D =n 2 + W 2 + [Jd^W + d^WJ] (5.12) 

In the non-relativistic limit, II = 0, and equation (|5.12p reduces to the 
quantum potential for the Pauli particle, pQ, 

Q P = W 2 + [S{VW) + (VW)S] (5.13) 

where 2pS = ipi^^^R is the non-relativistic spin limit of J. W is the 
non-relativistic limit of W^ 1 . 

The pseudoscalar part of equation (15,101) is simply [Jd^P^ — d^P^J] = 0. 
This puts a constraint on the relation between the spin and the momentum 
of the particle. In the non-relativistic limit this term vanishes. 

5.3 The Time development of the Spin 

Let us now turn our attention to equation (j5.4|) and show that it leads to 
an equation for the time development of the spin of the Dirac particle. By 
substituting equations (|5.8p and (|5.9p into equation (|5.4p we find 

J ■ -P-W + JA d^W^ = 

where we have written 

2j-d^ = Jd^ + d^j 

2P-W = PW + WP 
2 J A d^W^ = Jd^W» -d^W^J. 

All of these terms are Clifford bivectors so that equation (|5.4p gives just one 
equation. We can now simplify this equation since 

p {p . w) = -d» P {p^ ■ j) - p (p^ ■ d»j) 
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so that 



d^pP") ■ J + p(P^ ■ d»J) + P {J A d^W) = 

However since 2pP^ = T^°, the conservation of the energy- momentum ten- 
sor implies 

0„(r"°) = 2d^ P p») = o 

so that we have finally 

P M - d»J + J A 8^ = (5.14) 

This equation describes the quantum torque experienced by the spin of the 
particle in the absence of any external field. Coupling to an external field is 
achieved in the usual manner by replacing by — ieA^. The equation 
(|5.14p reduces to the quantum torque equation for the Pauli particle [lj 

(d t + — -)S = — (VW A S) 

m m 

Here P is the three-momentum and S and W have the same meanings as in 
equation (|5.13p . 

6 Conclusions 

In this paper we have presented a complete description of the Bohm model 
of the Dirac particle for the first time. This result demonstrates that the 
common perception that it is not possible to construct a fully relativistic 
version of the Bohm approach is incorrect. For simplicity we have considered 
only a neutral Dirac particle in order to bring out the novel quantum features 
of this approach. The extension to the case of a charged particle is however 
straight forward in principle and will not been presented here. 

In detail, we have obtained expressions for the Bohm energy-momentum 
density and for the quantum potential for the neutral Dirac particle. We 
have shown that these quantities become identical to those found for the 
Pauli particle first presented in [6] and recently re-derived in more general 
terms of the Clifford algebra £3,0 over the reals in [I]. We have also obtained 
the fully relativistic expression for the time evolution of the spin components 
of the Dirac particle. We have shown how this reduces to the non-relativistic 
version presented in both the above papers. The numerical details of this 
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non-relativistic version have been presented in Dewdney et al. [13] • We will 
examine numerical details of the relativistic case in another paper. 

It should be noted that the approach to the Dirac particle presented 
here is different from that presented in Bohm and Hiley [4| and in Doran 
and Lasenby [13], [33]. In these approaches it is the Dirac current 4>l1o4>r 
that is used to calculate trajectories. 

It should also be pointed out that Horton, Dewdney and Nesteruk [15] . 
although not motivated by the formulation of the theory in terms of Clifford 
Algebras, have proposed to base a Bohm- type theory on the use the eigen- 
vectors and eigenvalues of the energy momentum tensor to define particle 
trajectories and densities for massive relativistic scalar particle using the 
Klein-Gordon equation [15] and vector bosons [16] . They also give explicit 
examples of the different trajectories associated with the current (which 
yields space-like motions) and with the flows of energy momentum (which 
are necessarily time-like) for the case of scattering from square potentials. 
A similar approach has also been employed to develop a covariant version 
of the Bohm field ontology approach to quantised scalar fields [IT] . 

In this paper we have again emphasised that the Dirac current is different 
from the Bohm energy-momentum current. This implies that, in general, 
it possible to generate two different sets of trajectories. The appearance of 
two sets of trajectories is a novel feature of the relativistic domain. This 
difference is not unexpected as the earlier work of Takabayasi [9 J has already 
anticipated this feature. So too did Tucker [35] who used the Kahler-Atiyah 
algebra, a somewhat different approach to Clifford algebras, to show how 
these trajectories arise from Killing vectors. What our results show is that 
the difference disappears in the non-relativistic limit so the earlier work did 
not anticipate this possibility of doubling. We discuss the reason for this 
difference but feel that a more detailed investigation is necessary. We will 
leave this discussion to another paper. 
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